General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



C£- isaaii 


TRW REPORT NO. 30584-6002-RU-00 

N79-22049 

Dnclas 
25234 

EVALUATION OF THE DISCRETE 

VORTEX WAKE CROSS FLOW MODEL 
USING VECTOR COMPUTERS 

Part H - Users Manual for DIVORCE 


(NASA-CK-152271) EVAIOATION OF THE DISCBETE 

VORTEX HAKE CROSS FLOW MODEL USING VECTOR 

COMPUTERS. PART 2: USER'S MANUAL FOR 

DIVORCE (TRH Defense and Space SystecDs 

Group) 70 p HC A04/MF A01 CSCL 01A G3/02 


F. D, Deffenbaugh and J. F, Vitz 


Prepared for 

Ames Research Center 

Under 

CONTRACT NAS2-9579 
JANUARY 1979 



TRW REPORT NO. 30584-6002- RU-00 


EVALUATION OF THE DISCRETE 
VORTEX WAKE CROSS FLOW MODEL 
USING VECTOR COMPUTERS 

Part II - Users Manual for DIVORCE 


F. D. Deffenbaugh and J. F. Vitz 
TRW Defense and Space Systems Group 
One Space Park 

Redondo Beach, California 90278 


Prepared for 

NASA/Ames Research Center 
Moffet Field, California 94305 
and 

U.S. Army Missile Research and Development Command 

Redstone Arsenal, Alabama 35809 

Under 

Contract NAS2-9579 
January 1979 



oereNSt am space systems cpoup 



TABLE OF CONTENTS 


INTRODUCTION 

1.0 PROBLEM DESCRIPTION . . . 

2.0 INPUT 

2.1 Card Input 

2.2 Input Descriptions . 

3.0 OUTPUT 

3. 1 Printed Output . . . 

3.2 Data File Output . . 

4.0 SAMPLE CASES 

5.0 OPERATION 

6.0 REFERENCES 


. 1 
. 2 
. 3 
. 3 

. 3 
. 13 
. 13 
. 14 
. 31 
. 67 
68 


1 



INTRODUCTION 


Volume 1 of this two part report presents the theoretical 
development of the discrete vortex wake cross flow method for 
predictinQ forces on bodies of revolution at high anQls of attack. Results 
for several missile test cases are also presented in Volume 1. This 
volume of the report is the users manual for the Discrete VORtex Cross 
flow Evaluator (DIVORCE) computer program. DIVORCE was developed 
in FORTRAN IV for the CDC 6600 and CDC 7600 machines. Opti.ual calls 
to a NASA vector subroutine package are provided for use with the 
CDC 7600. 
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1.0 PROBLEM DESCRIPTION 


The Discrete VORtex Cross flow Evaluator Program (DIVORCE) 
is used to predict the aerodynamic force and moment coefficients on bodies 
of revolution at high angles of attack. 

The time development of the flow behind an impulsively started 
circular cylinder whose radius is changing in time is computed using a 
large number of discrete point vortices to describe the wake region. 

The point vortices are introduced into the wake from boundary layer and 
rear shear layer separation points. 

The unsteady integral momentum boundary layer equations are solved to 
detennine boundary layer separation for early times after impulsive 
start and for times where the cylinder radius is changing. For times 
later in the solution where the cylinder radius is constant and quasi- 
steady flow may be assumed Stratford's (1) criterion is used to determine 
separation. A form of Stratford’s criterion is included for turbulent 
as well as the laminar flow case. The unsteady integral momentum 
solution is only valid for laminar flow. 

The viscous cross flow analogy (Ref. 2) is invoked to relate the 
two dimensional results to the distribution of normal and side force on 
the three dimensional body of revolution. The rate of change of the 
cylinder radius depends on the three dimensional geometry. 

The normal and side force distributions are integrated to yield the 
total normal and side forces, and pitching and yawing moments. 
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2.0 INPUT 


2.1 CARD INPUT 

Except for the case header card, all card input is accomplished 
via NAMELIST $ENGAGE. The case header card must be the first card in 
the input data deck. Up to 60 characters may be input starting in 
column one. Following this card is the $ENGA6E card and then all of 
the input items for this case. The input is terminated by a $END 
card. 

2.2 INPUT DESCRIPTIONS 
2.2.1 Body Geometry 

DIVORCE is able to operate on any axisymmetric body. A standard 
ogive cylinder is the nominal body shape and the program input 
needed to describe it is minimal. Any other body must have radius given 
as a function of length. 

For an ogive cylinder the length and maximum diameter must be 
entered in dimensional values. The nose length may be entered as a 
dimensional value or as a number of diameters. If the second option 
is chosen, the nose length is multiplied by the diameter prior to 
the calculations. Given these values a tabl<^ of radius versus length 
is generated. 

For a general body this table of radius versus length must be 
entered. There are a maximum of 51 values in each table. 

Given these tables, tables of time, nondimensionalized length, 
nondimensionalized radius and its first and second derivatives, and 
values of the equivalent two dimensional cylinder radius, a|^, and its 
first and second derivatives are generated. Some of tnese parameters are 
used in the final calculations for the aerodynamic coefficients. 

For the standard ogive cylinder, the values of a|^, a'j^ are 
computed at each time step. For other bodies, these values are 
found by linear interf )lation upon the created tables. 
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2.2.2 Execution Controls 


An initial distance along the body length is given at which to 
start execution. A value of time is found for this distance. A' 
time step is given either as a constant or variable. If DELT is 
entered it is assumed that this is a constant time step. For a 
variable time step a minimum and maximum step size are given and 
DELT will vary between them. DELT must not be entered to use this 
option. Nominal values for the limits are .05 and .125, The 
program will continue until one of the following conditions has 
been fulfilled. 

1) The number of time steps = KFINAL 

2) The length along the body ^ ZFINAL 

3) The amount of CPU time ^ TFINAL 

2.2.3 Restart Capability 

It is possible to stop the execution at an intermediate point, 
check the answers, and continue if desired. If RITE0UT = 1 , the 
common blocks will be written into a file (TAPES) when one of the 
three conditions in Section 2.2.2 has occurred. This file may be 
stored on disk or tape and will be read in at the beginning of the 
next run. READIN must be set to 1. and the data will be read from a 
file (TAPE4). The restart capability can be used any number of 
times, but care must be taken to update the card input data for 
KFINAL, ZFINAL and/or TFINAL. Most of the input data is not saved 
and must be reentered on the restart run. 

2.2.4 Perturbation 

For a certain distance along the body it may be desired to per- 
turb the value of the vortex strength parameter, (SIGMA). ZPERT and 
ZPEND define this segment. SIGMA is perturbed by an amount which 
is input as a fraction of SIGMA. For example if DELSIG = .1 the 
change in SIGMA is defined by 
Aa = (-.1 )a 
and 

o' = 0 + Aa in the range ZPERT to ZPEND 
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2-2.5 Pressure Calculatinng 


The pressure calculations are done every KPRES time steps. 

Values of and and an array of pressure on the body as a 
function of the angle around the body are computed. These may 
be output as files or listings as desired. A time history of 
''d stored for use in the computation of CN, CY, Cf„ and 

C|^ at the end of the run. ZLAMDA is used in the computation of 
^m *"n nominally zero. 

2.2.6 Coalescence 

It is possible to combine two or more vortices into one if 
they are sufficiently close together. Nominally no coalescence 
will occur but if KRC0AL and/or KC0AL are set to positive integers, 
coalescence will occur every KRC0AL or KC0AL time steps. A value of EPS 
can be entered (nominally = . 1 ) to use in the detemiination of the 
proximity equations. If two vortices satisfy the proximity test 
they are conbined to form one. This will cause differences of some 
degree in the answer, but will decrease execution time significantly. 

2.2.7 Vectorization 

Both loop code and vector code has been used in the program. 

At a low number of time step the loop code executes faster than the 
vector code. But as the number of steps nears 40 the execution is 
about the same and then the vector code becomes faster. An input 
parameter (KVEC, nominally = 50) is available which tells the program 
to shift to the vector code when the number of time steps reaches it. 
There appears to be a significant round off difference in the vector 
routines which will cause somewhat different answers. 

2.2.3 Boundary Layer Separation 

Two different solutions to the boundary layer equations are 
used to determine the separation point. For times less than t = 4.0 
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the unsteady integral momentum equations are solved for laminar 
flow. For times greater than t=4.0 , and provided that the 
cylinder radius is constant Stratford's method Is used to 
determine the boundary layer separation points. Stratford's method 
may be used for either laminar or turbulent flow. If LAM = 0 the 
laminar scheme is used. If LAM = 1 Stratford's turbulent equations 
are used and the transition length (XTR) must be entered. 

2.2.9 Print Flags 

Various types of printed output exist in the program. The 
minimal amount of print will normally be used, especially for long 
runs. Each time step consists of one line of print (see Section 3.0). 
The vortex information may be printed at every KPRINT steps and the 
pressure data may be printed every KPRINTP times that the calcula- 
tibns are made. In addition, intermediate information in a number 
of areas of the program may be printed for use in debugging. Care 
should be taken in using these print flags because a lot of print 
can be generated in a very short time. 

2.2.10 Plot Flags 

It may be desired to output data at specific time steps onto a 
disk file or tape which may then be run through a plotting program. 
Three distinct types of output exist and their formats are described 
in Section 3,2. Briefly they are the vortex information, pressure 
information and aerodynamic coefficient information. 

The vortex information may be written out at every KL0UT time 
steps or a table of nondimensional ized body length values [Z0TV (100)] 
may be input. If one method is used the other can not be used. For 
the second method, the data is written at the first time step past a 

value in the table. The output will terminate upon encountering a 
zero in the table. 

The pressur^' information is written out at the first time step 
for which it is calculated past a value in the table Z0TP(1OO). 

Likewise, output will terminate upon hitting a zero in the table. 
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The aerodynamic coefficients are written on the same file as the 
pressure terms for any run with KPRES not equal to zero. It is written 
regardless of the existence of any pressure data. 
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TABLE 2.1 


SYMBOL 

ISYM 

LAM 

RITE0UT 
READ IN 


$ENGAGE NAMELIST INPUT 


DESCRIPTION 


* * * FLAGS * * * 

Symmetry flag for the pressure 
calculations. 

= 0 asyiPTetric pressure integration 
= 1 symmetric pressure integration 

Boundary layer separation 
flow flag 

= 0 for laminar flow 
= 1 for turbulent flow 

Restart flag 
= 0 for no restart 
= 1 to write restart data on file 
TAPES for use with a continu- 
ation run. 

Restart flag 
= 0 for no restart 
= 1 to read restart data from 
file TAPE4 in order to begin 
execution at last point of a 
prior run. 


NOMINAL 

0 . 

0 . 

0 . 

0 . 
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SYMBOL 

DESCRIPTION 

NOMINAL 


* * * CONSTANTS * * * 


RC 

Discrete point vortex core radius 

0.05 

XTR 

Transition length in radians from 0 = 0.0. 
Used only if LAM =1 

0.0 

ATTACK 

Angle of attack in degrees 

0.0 

SIGMA 

Vortex strength parameter 

0.6 

DELSIG 

Amount to perturb SIGMA in the range 
ZPERT to ZPEND. Entered as a fraction 
of SIGMA. DELSIG (comp) = -DELSIG (input) 
*SIGMA and SIGMA (pert) = SIGMA (input) + 
DELSIG (comp) 

0.1 

ZLAMDA 

Moment center measured from the nose 
and nondimens ionali zed by body length. 
(A = Z*/£). 

0.0 

EPS 

Coalescence radius 

0.1 


* * * BODY GEOMETRY * * * 


OGIVE 

Ogive cylinder flag 
= 0 for general body shape 
= 1 for ogive cylinder 

1.0 


* * OGIVE CYLINDER * * 


ZLENGTH 

Total body length 

0.0 

DEE 

Maximum body diameter 

0.0 

ZLN 

Nose length 

0.0 

ZLND 

Nose length as a number of diameters 
ZLN = ZLND * DEE 

0.0 
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SYMBOL 


DESCRIPTION 


NOMINAL 


ZSTAR 

RSTAR 

KFINAL 

KPRES 

KVEC 

KC0AL 

KRC0AL 

TFINAL 

DELT 

DELTMIN 

DELTMAX 


* * GENERAL BODY * * 


(51) Argument table of body length values. 

The last value is set to ZLENGTH. 

(51) Function table of body radius values. 

The largest value is used to compute 


* * * TIME STEP CONTROLS * * * 


The number of time steps at which 
to terminate execution. 

Pressure calculations are made at 
every KPRES time steps. If KPRES = 0 
no calculations are made. 

Vector method flag. When the step 
counter reaches KVECj the vector routines 
are invoked. 

Coalescence of all vortices occurs everv 
KC0AL time steps. 

Coalescence of rear vortices only occurs 
every KRC0AL time steps. 


* * * TIME CONTROL * * * 


Execution is terminated when the CPU 
time reaches TFINAL. 

Constant value of delta time for each 
step. 

Limits of delta time as it varies throughout 

the run. DELT must not be entered to use this 
option. 


0 . 

0 . 

512. 

5. 

50. 

0 . 

0 . 

500. 

0 . 

.05 

.125 
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SYMBOL 

KPRINT 

KPRINTP 

KL0UT 
Z0TV (100) 

Z0TP (100) 

ZFINAL 

ZINIT 

zpert) 

ZPEND( 


DESCRIPTION 


* * * OUTPUT FLAGS * * * 

Vortex information is printed 
at every KPRINT time steps. 

Pressure values are printed every 
KPRINTP points, assuming that they 
have been calculated. That is 
whenever 

M0D (K,KPRES) = 0 and 
M0D (K, KPRINTP) = 0 

Vortex information is written onto 
an output file (TAPE6) at every 
KL0UT points. If Z0TV is being 
used this must not be entered. 

Array of nondimensionalized body 
length values at which to write 
vortex information on the output file 
(TAPE 6). The data at a point immediately 
after a value of Z0TV is written. 

Array of nondimensionalized body 
length values at which to write the 
pressure values around the body on 
an output file (TAPE7). The first 
point at which the pressure values 
have been calculated past a value of 
Z0TP is written. 


* * * BODY LENGTH CONTROLS * * * 

Nondimensionalized value of the 
length on the body at which to 
terminate execution. 

Nondimensionalized value of the 
length on the body at which to 
initialize execution. 

Nondimensionalized values along 
the body which define the beginning 
and end of the perturbations on SIGMA. 


NOMINAL 

0 . 

0 . 

0 . 

0 . 

0 . 


1 . 

0 . 

0 . 

0 . 
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Sli'IBpL 

DESCRIPTION 

NOMINAL 


* * * DEBUG PRINT FLAGS * * * 


PRV 

Nonzero to print rear vortex 
separation information. 

0. 

PBL 

Nonzero to print steady boundary 
layer separation information. 

0. 

PAP 

Nonzero to print unsteady boundary 
layer separation information. 

0. 

PC0AL 

Nonzero to print coalescence 
information. 

0. 

PFIX 

Nonzero to print vortex motion 
corrections and pressure calcula- 
tion corrections. 

0. 
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3,0 OUTPUT 


3.1 PRINTED OUTPUT 

3.1.1 Nominal 

For each execution of the program, the following items are 
written to the print file. Page one contains all the NAMELIST 
input parameters and their values for the case except for the 
arrays Z0TP, Z0TV, RSTAR and ZSTAR and the body geometry data, 
which appears on page two. Listed therein are the input values 

and the calculated values. The definitions of these values are 
given in Table 3-1. 

Following the input and geometry information the program 
lists one line for each time step. These values are defined in 
Table 3-2. The last page contains a summary of the aerodynamic 
coefficients and the separation angles. These values are defined 
in Table 3-3. The tabular values are stored at each pressure 
calculation step and the scalars are computed at the end of the 
run. The local normal force coefficient and local side force 
coefficients are obtained from the two dimensional drag and lift 
values in Table 3-2 by the factor shown below. 

= CDPK • sin2« (2AW)/DEE 

Samples of the print format for these pages are given in 
Figures 3-1, 3-2, 3-3 and 3-4, 

3.1-2 Vortex Information 

At each KPRINT time point the information for all of the 
vortices is printed. The front top and bottom vortices are printed 
first followed by the rear top and bottom vortices. The parameters 
are described in Table 3-4 and the format shown in Figure 3-5. 

3.1.3 Pressure Information 

The pressure information and drag and lift coefficients are 
computed every KPRES times. These values are printed whenever 
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M0D (K, KPRINTP) = 0 and M0D (K, KPRES) = 0. The parameters are 
defined in Table 3-5 and the format given in Figure 3-6. 

3.1.4 Coalescence Information 

Whenever the vortices are coalesced, information concerning 
this function can be printed out. The vortex information is printed 
both before and after the coalescence and a list of the actions 
taken is given. The longer of the two tables is always listed on 
the left side. The sets of vortices are referred to by the following 
numbers 

1 - Front top 

2 - Front bottom 

3 - Rear top 

4 - Rear bottom 

Table 3-6 and Figure 3-7 define and describe this information. 

3.1.5 Debug Information 

Information concerning the rear vortex separation, the unsteady 
boundary layer separation, quasi -steady boundary layer separation, 
and some fixup information can be printed on user option. Care 
should be exercised so that a ridiculous amount of print is not 
generated. The parameters printed in these sections of the program 
can be traced through the subroutine listings. 

3.2 DATA FILE OUTPUT 

3.2.1 Vortex Information 

The vortex information is written onto a data file either at 
every KL0UT time points or at a specified value of body length given 
in table Z0TV. This data is written on TAPE6 and has the following 
format: 

One logical record consisting of N words where 

N = 4 + 4K0 • 6 

K0 = MAX0(KT, KB, KRT, KRB) 

KT = number of top vortices 
KB = number of bottom vortices 
KRT = number of top rear vortices 
KRB = number of bottom rear vortices 
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The data is the time step number (K), the time (T), the length 
along the body (Z) and K0 followed by K0 sets of a counter (I), 

X and y position, x and y velocity and vortex strength gamma for 
each of the four types of vortices. 

The F0RTRAN write statement used is: 

WRITE (6) K, T, Z, K0, (I,(X(I,J), Y(I,J), XD0T(I,J), 

YD0T(I,J), GAMMA J = 1, 4), I = 1. K0) 

3.2.2 Pressure Information 

At the first point at which the pressure information is 
calculated past a value in ZOTP, the equivalent three dimensional 
pressure as a function of angle around the body is written on a file 
(TAPE7). Each logical record contains 126 words which are the number 
of words to follow (NITS = 125), the step number (K), time (T), the length 
along the body (Z) and 61 pairs of body angle (PHI) and pressure (PRES). 
The F0RTRAN write statement is: 

WRITE (7) NITS, K, T, Z, (PHI(I), PRES(I), I = 1, 61) 

At the end of the run a dummy record of the same length but with 
K set to zero is written. This is used to signify that the next record 
contains the drag and lift coefficient data. 

3.2.3 Drag and Lift Data 

If any pressure information was calculated the last record on 
file TAPE7 will contain the drag and lift coefficients as functions 
of body length. This record will appear even if no output was generated 
in Section 3.2.2. The logical record is of length N where N = 3 ’ NPC + 1 
where NPC = number of times the pressure calculations were made. 

The record consists of the number of words to follow (NITS = 3 ' NPC) 
plus NPC sets of body length, drag coefficient and lift coefficient. 

The F0RTRAN write statement is: 

WRITE (7) NITS, (ZZ(I), CSD(I), CSL(I), 1=1, NPC) 
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TABLE 3-1 

BODY GEOMETRY DATA 


SYMBOL 

LENGTH 

DIAMETER 

AREA 

F 

AW 


ZSTAR 


RSTAR 

Z 

T 

RZR0 

DRZR0 

DDRZR0 

A 

AD0T 

ADD0T 


DESCRIPTION 


Overall body length as input as ZLENGTH (OGIVE = 1) 
or as extracted from table ZSTAR (0GIVE = 0.), £ 

Maximum body diameter as input as DEE (0GIVE = 1) 
or as extracted from table RSTAR (0GIVE = 0.), d 

Frontal area, S=ird2/4 
Fineness ratio, z/d 


Characteristic body radius. 

non-maximum diameter 

ended bodies, ie., ellipsoids, boatails, etc., AW 
IS computed as the average of the body radius, a 

Table of lengths along the body. For 0GIVE = 1 the 
value of ZLENGTH is divided by 50 to obtain the’ 

injuf ° 


Table of radius values along the body. For 0GIVE = 1 
these are computed. For 0GIVE = 0. this table is input. 


Nondimens ionalized values of body length, Z*/z 
we'body! to the length along 

Nondimensionalized values of body radius, R = 2R*/d 


First derivative of radius with respect to length. 

Second derivative of radius with respect to length. 

Circular cylinder radius at time (T), a(t) = - R fZl 
For 0GIVE = 1, A = RZR0 ^ 2a 


First derivative of A with respect to time. 
Second derivative of A with respect to time. 
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TABLE 3-2 

CIRCULAR CYLINDER TIME DEPENDENT RESULTS 


SYMBOL 

DESCRIPTION 

K 

Number of current time step. 

TIME 

Time, t = y|- 

Z 

Nondimensionalized value of body length, Z 

THETAS 

Top separation angle in degrees. 

THETASB 

Bottom separation angle in degrees. 

THTASR 

Top rear separation angle in degrees. 

THTASRB 

Bottom rear separation angle in degrees. 

CDPK 

Drag coefficient. 

CLPK 

Lift coefficient. 

AK 

Cylinder radius A at time step K. 

AKD0T 

First derivative of A at time step K. 

AKDDT 

Second derivative of A at time step K. 


= 
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TABLE 3.3 


AERODYNAMIC COEFFICIENTS 
SYMBOL DESCRIPTION 


1 

Nondimensionalized body length, l*/s. 

CD 

0-, 2 

Local normal force coefficient, ^ Sin a CDPK 

CL 

Local side force coefficient, ^ Sin a CLPK = C 

THETAS 

Top separation angle in degrees 

THETASB 

Bottom separation angle in degrees 

CN 

Normal force coefficient, N/q S 

^co 

CY 

Side force coefficient, Y/q S 

CSM 

Pitch moment, M/q Sz 

*co 

CSN 

Yaw moment, MY/q^Sz 
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SYMBOL 

N 

K 


X 
Y 

XD0T 
YD0T 

XDDT 
YDDT 

GAMMA 

"B" suffix 

"RT" suffix 

"RB" suffix 


TABLE 3-4 
VORTEX INFORMATION 

DESCRIPTION 

Point vortex number, x^, y^, t^ 
Current time step counter. 

Top vortex position. 

Top vortex velocity. 

Top vortex acceleration. 

Top vortex strength. 

Designates bottom vortex data. 
Designates rear top vortex data. 
Designates rear bottom vortex data. 
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TABLE 3-5 


SYMBOL 

T 

Z 

AK 

AKOOT 

THETAS 

DEG 

PHIVT 

PHIPT 

2(PHIT) 

-PSKSQD 

CPK 

CP3D 

UTAN 


PRESSURE INFORMATION 


DESCRIPTION 


Time. 

Distance along the body. 

Cylinder radius at timestep K. 

Derivative of cylinder radius at timestep K. 
Separation angle. 

Angle around the body. 

Component of the potential term in the Bernoulli 
equation due to point vortices, 

Component of the potential term in the Bernoulli 
equation due to changing radius, 

Time derivative of the potential term in the 
Bernoulli equation, + 4 >P,^) 

Negative of the square of the surface velocity. 

Two dimensional circular cylinder coefficient of 
pressure at timestep K, (2$,. +1 -u^) 

V 0 

Three dimensional surface pressure, {2a/d)*Sin^a*CPK 
Surface velocity, 
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TABLE 3-6 


COALESCENCE INFORMATION 


SYMBOL 


N 



"B" suffix 

11 

12 

LISTl 
LIST2 


DESCRIPTION 

Vortex number, x^, y^, 

Position and strength of vortex of first 
set of data. 

Designates second set of data. 

Locations of vortices being combined. 

Identification of type of vortices 
being combined. 
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UuCU(<tNT SAHPLt CASE 


K^•1^AL 

KPMES 

KVPC 

KCQAL 

KRCUAL 

KLOUT 

ttPiUfiT 

KPRlfJTP 

isvn 

LAM 

Ao 

5 

^r\ 

0 

C 

c 

C 

0 

u 

0 

Rt ADIN 

RITEOUT 

DQV 

PUL 

PAP 

Pf IX 

PCOAL 

OGIVE 



0 . 

0 . 


0 . 

0 . 

0 . 

0 , 

1. 



Tf INAl 

CELT 

O'^LTHIN 

OFl fttAX 

4FiMAL 

ZiNlT 

ZPLKT 

ZPENO 

RC 

XTR 

dO J.uOu 

0,000 

, 0^0 

,1<!5 

X, UbC 

.010 

.050 

, d 3 J 

• 05o 

0.000 

AATACK 

SI6«A 

0P| $TG 

7LAnDA 

cPS 






30 « obo 

.600 

,100 

0,060 

.100 







FIGURE 3-1 NAMELIST INPUT PARAMETERS 
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4.0 SAMPLE CASES 


Three sample cases are Included. The first is a nominal run 
with 40 time steps. The second is the same run with coalescence 
at time step 35. The third contains the vortex and pressure 
information printouts. 
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.267722 

-lfc0,680 

1.759928 

.081762 

*373610 

.236327 

-160.360 

1. 7599 26 

.081762 

.897896 

.210366 

-lfc0,600 

1,682651 

•082228 

.920127 

•Ibbalb 

-160,506 

1.682651 

•082228 

.960059 

,la0799 

-157,207 

1.682651 

.082228 

.937658 

. 135360 

-162,337 

1.682651 

.032228 

.972108 

.109699 

-160.161 

1.682651 

.032228 

•933326 

.083337 

-168.123 

.866121 

.181671 

.992667 

.058918 

-160.151 

.865121 

•181671 

,997857 

•U303U7 

-162.365 

.866121 

.181671 

,999970 

.003572 

-162.663 

•865121 

• 1816 71 

1.000000 

O.OOOOOO 

■“158.667 

*865121 

•191671 

1.000000 

O.OOOOOO 

-167.360 

■562266 

-.010532 

l.OOODVO 

0.000000 


AKUDT 

".2607t)6 

-•258{Jb2 

-,2!>ru3!» 

-.255ZA1 

-.2535QA 

-.251816 

-.250iaU 

-.2A859Z 

-.ZA7U5Z 

-,ZA555f 

-.ZAAllZ 

-.2VZ591 

-.ZA1ZZ7 

-.Z197ZS 

■>.231191 

-.236030 

X.235UA9 

-.233557 

-.Z31D61 

-.230Z7Z 

-,ZZ869f 

-.ZZTlJZ 

-.225053 

-V22518J 

-.222783 

-.2Zl5:iA 

-.220208 

-.219056 

-.218001 

“.ziruM 

-.216263 

-.215586 

-.215062 

-.216065 

-.2143V6 

-.216299 

O.UbOOOO 

U.OOOOOQ 

o.oobouD 


T 

CD 

CL 

.C26P13 

.70o23<? 

O.OOOOOC’ 

.047*^40 

.3ng531 

O.OOUOOO 

.069459 

.391502 

.003236 

,09f534 

.476160 

.012845 

.135297 

.439962 

.U20436 

.179761 

.370613 

.020457 

.229937 

. 216260 

.045410 

.2P2390 

.136061 

-.002645 

CN 

C7 

CSM 

l.OFf-350 

.056474 

-.167943 


fH(:TAS 

112. 260 
10B.566'«d2 
Ufa, 207993 

i01.7l!j370 
98.6ii8 7ub 
97.^37331 
92.203964 

CSN 

-.011J62 


rrffcTASD 

-112.260719 

-108.536482 

-106.441760 

,7105.015390 

-102.573305 

-100.719952 

-90.464343 

-93.956950 



(namelist input) 


JVTXIN 

KFINAL*40 

ZPERT*« 05 

ZPEND=,Z5 

AATACK-30 

ZLENGTH«78.272 

ZLND-2.598 

D£t*7«6 

ZIMT-.Ol 

KVEC*30 

KC0AL»35 

KRCQAL-35 

PCOAL-1 

SEND 


CASE 2 


38 


DOCUMENT SAMPLE CASE 


KFINAL 

KPRES 

KVEC 

KCOAL 

KRCOAL 


5 

30 

3£i 

35 

RtAUiN 

RITHOUT 

PRV 

PBL 

PAP 

0. 

0. 

0. 

0. 

0, 

TFIHAL 

delt 

DELTHIN 

OFLTMAX 

zfinal 

S'itJmOQO 

QtOOO 

• 0?Q 

• 125 

1.000 

AATACK 

SIGMA 

DELSIG 

ZLAHDA 

EPS 


KLOUT 

KPAINT 

KPRlMTP 

ISTH 

LAN 

0 

Q 


0 

0 

PFIX 

0. 

pcdal 

1. 

□GIVE 

1. 



nun 

IPcRT 

ZPeNO 

. RC 

XTK 

• 010 

.050 

<250 

.050 

o»ooo 



o 


LENCTH 

76,2720 


iilAk pstar 


O.OUOO 

.0000 


.5i973 

3.1309 

1.1391 

4 .0963 

1.6272 

6. ..615 

2.0630 

7,0272 

2.4479 

9.392D 

2.7830 

10. 95 6 j. 

3.0693 

12.523t 

3.3076 

14.0U90 

3,4965 

15.0044 

3.6425 

17.2l9t 

3.7400 

16. 705 3 

3,7913 

20. 3507 

3.8 000 

21.9162 

3.6000 

23,4blo 

3.8000 

25.0470 

3.0000 

20.6125 

3.0000 

23.1779 

3.8000 

Z9./434 

3.8000 

si.juau 

3.0000 

32.6742 

3.8000 

34,4397 

3.HOOO 

3o. Ou5 1 

3.8000 

37,5700 

3.8000 

39 . 1 3ou 

3.0000 

40.7014 

3.8000 

42.2669 

3.8000 

43.8323 

3,0000 

45,3978 

3.8000 

46.9o32 

3.8000 

48.5266 

3.8000 


nUHtTfcR 


7«6UU0 
7 T 


0.0000 

O.UOOU 

.0200 

■ 2378 

*0400 

.4757 

.0600 

.7135 

.0800 

.9514 

.1000 

1.1692 

.1200 

1.4271 

.1400 

1.6649 

.1600 

1.9020 

• 1800 

2.1406 

.2000 

2.3764 

.2200 

2.6163 

.2400 

. 2.6541 

.2600 

3.0920 

.2000 

3.3298 

.3000 

3.5677 

.3200 

3.6055 

.3400 

4.U433 

.3600 

4.2812 

.3800 

4.5190 

.4000 . 

4.7569 

.4200 

4.9947 

. 4400 

5.2326 

.4600 

5 »4 7q4 

.4800 

5, 7003 


• ^DilO 9 * 9 't 6 i 

•5200 6,i839 

•5^00 fi t 4216 

•5600 &*u696 

•5500 fi.a976 

,6000 7tl3i>3 

•6200 7,3732 


DbCUHfiMT SANRL£ 
WD7 GEOMcTRT 
AR£A 


45 (3646 
K 7 K£i DRZRO 


.0000 

6.2333 

.1572 

7.4900 

• 2998 

6,7716 

• 4282 

6.0749 

• 5429 

5.3969 

• 6442 

4,7348 

.7324 

4*0665 

.8077 

3.4496 

• 6704 

2.8222 

• 9207 

2.2024 

.9566 

1.5065 

.9842 

.9788 

.9977 

.3716 

I.UUOO 

U.OOOO 

1,0000 

O.OuOO 

l.buOO 

0.0000 

1, uooo 

0.0000 

l.UOOO 

0.0000 

1,0000 

0.0000 

IfOOuO 

0.0000 

1,0000 

D.OCOO 

l.OOOC 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

1.0000 

0.0000 

l.uOOO 

0. 0000 

1.0000 

0.0000 

1.0000 

0,0000 

l.uooo 

O.oooo 

1.0000 

0.0000 

i.OOOO 

0.0000 

1.0000 

0.0000 


CASE 


F 

X0*29A9 

UORZRQ A 

-37.5531 #0000 

-36.5l2^ .1572 

-35.3507 .2998 

-34«34£>ip .4282 

-33.4007 ,5429 

-32»740fi »6442 

-32.113P ,7324 


-31,6909 

.8u77 

-3l,ib44 

,6 704 

-30,8262 

.9207 

-30.6777 

«9?{>o 

-30.4097 

.9842 

-30.3217 

.9977 

o.ooou 

l.UUJU 

O.OUOO 

l.OuUO 

O.UuOO 

i.u Jou 

O.OUOO 

1.0 uuu 

0. uuOO 

l.C JuJ 

u.oooo 

1 .OUilO 

O.OUOO 

l.U jOj 

O.OUOO 

l.uooo 

u.oooo 

l.OUljo 

O.OuOO 

I.OOOO 

0.0000 

I.OOOO 

0.0000 

i.uouu 

O.OouO 

l.OUUO 

O.vOOO 

l.UUu J 

o.oooo 

I.UUOO 

0.0000 

1. 0 ju '1 

O.OOUO 

l.U JO 1 

o.oooo 

1 . 0 Ju J 

O.OUOO 

1.0 uu ) 



' AOOT 

*6923 

• 6 296 

• 5694 
.5108 
»4£36 
«39ol 
.3436 

• 2901 
.2373 
.1652 

• x336 

• 6333 
.0312 

d.O Jou 
JfuOliO 

O.COLO 

ViCusjO 

J.'jOOO 

OiOOuO 

O.uOCiO 

\i>uCOO 

0.0000 

u.OOUO 

w .OOOu 

Otoool 

OlOuuo 

OtObOU 

U* OOuO 

0.0 ooi) 

1 • JOOO 
5 • L tlOu 

O.vOOO 


AH 

3.8000 

AOOUT 

-.2677 
-.2502 ' 

-.2500 
-.2429 
—•2367 
-.2315 
-.2371 
-.2234 i 
-.2204 I 
-.2180 ; 
-.2102 
-.2100 : 
-.2144 
v/.OubO 
O.OuOO 

o.oooo 
0.0000 ' i 

O.OOCO ; 
O.CbOU 
O.UoOO 
0. JbOO 

o.cclou j 

O.OCiOU 

0.0000 I 

O.uuUCi 
0.0000 
O.OOUO ; 
O.OgoO 
O.OuCO 
U.UUOO 
O.OuOO 
0.0000 


‘Sociaaoaaaooaao 
Aoosooooooooso 
50533003333^00 
o o o .0 o o o o o o a o o o 

• «•»••• • m « • • •« 

00 '^Q"» 00 ' 3000 aO '3 


^ cao a ^ 

3 3 3 0^3 
3 0 0 0 0 

• • * « * 

O 00 41 3 


03 * 3^00000 

000000300 

OOOOOOOOO 

o aooooooo 

*•««««•«• 

ooooooooo 


9 -O 

9 30 

.0 

9 9 

9 0 

0.9 0 

•3 

9 9 

3 9 

3 0 3 

9 

0 0 

0 0 

0 0 9 

9 

9 3 

• m 

• • • 

» 

• • 

9 •> 

3 0 3 

0 

0 9 


o o 
o o 
o 9 

oo 

# • 

o o 


OOOOOOOOOOOOOOOOOOO 
22 94 > 0 Q 0 0 0 00 90 309009 
09900 0 900 90 09 390 390 
99 990 0099 90 OOO 3 9 900 
*••••••« ••o. •••• 


«i««ipww»wwow 900909<9090 
O^OOrO’O'OOO 3 30 9000 3 09 
9 0 9 0 0 90 9 9 9 9 0 0 9 0 0 9 9 9 
090 OOOOo 09090. 3 0 0 30 3 

O 00 000 990 3 00 90 900 9 9 


00000000009 
00000909900 
00 0 900 90900 
00090 oooooo 


90000000 

OOOOOOOO 

00009909 

oooooooo 


0 9 0 0 0 9 0 0 9 O O O O 


9 O 9 O O O 


0000 00 9 0.0 O O 
O *9 099^000 30 
0 9 9 0 0.9 0 0 9 3 e 
9 9 0 9 9.0 9 9 9 3 9 



o o a o o o 

3 0 0 9 C3 O 
9 93 939 

0 0 9 0-30 



o o 

9 O 
0 O 
9 9 


f«4 9 0 f\4 o 9*TVfri -HJ^r^ftcno-oo^rtvi 

^ •r t> N 9 o fT> U 3 « C\J ^ o 4 ^ rs. . 3 * 

45 3 O tfl ^-50 a 


^ « O’ O*. O CT^ ^ e o O O ^ ^ ^ 


OCCOO^OCCO OC^COCQCCei 

ooooooccco ec^c o<-ooco 
^ 4 ^VC^|^r« 0 cPO(v*^ 4 ; 7 ■O(^‘^r’'Cc :9 
9 Kf’-r- tr a^rp cc cr o O’ cyfr C^O 


OOOOOOOOOOQOOOOOOOO 
OOOOOOCOOOO 000 30000 
o o o o o o a o o 0 3 o o o o o a o o 
csjctnooccsjct^ccccc oocb^cDaQa^viu 

• ♦•••••» « « • « • •« •« 

(nfnfvtcnrncnntPtr) cn »t> rn «ti <n tn m 4*1 


“"♦ O •!* 3 > *n P- OU 0 ,*> 1 > 'T 9 .M ^- -.H 9 

spgh^O- 4 %*^ 9 <>Jf*^cri 333 * 4 ' 9 ’ 4 \*-* o»V 

T ^ 

99 Ai.Nim> 4 * 3 « 0 -<*’"-’» a*r 

* « • •' .« • «r • • B « • "*' 

O 4> P% t3» ^ I'j 'Q 9 (M 41 -0 

>>A‘A 4 \\A 4 \ 4 V-a 0-0 
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K 

riMe 

Z 

TH^TA^ 


1 

.119 

.010000 

160,000 

- 130 . OOU 

2 

.169 

.014204 

140.399 

- 149,399 

3 

.219 

.018409 

117,417 

- 117.^17 

4 

.269 

*022613 

114.720 

— 114.720 

5 

.319 

■026818 

112.261 

- 112 , 26 i 

6 

.369 

,031022 

111.530 

- 111.530 

7 

.419 

.035227 

110. 786 

- 110,786 

8 

.469 

*039431 

110,147 

- 110 , j .47 

9 

*519 

.043635 

109,259 

- 109.269 

lO 

.569 

.047840 

104,586 

— lOO.iSt 

11 

.619 

.052044 

103.123 

- 103,123 

12 

.669 

.056249 

107.465 

- 107.465 

13 

.720 

.060508 

106.967 

- 107,053 

14 

.773 

.065041 

106.459 

- 106.6 25 

15 

.631 

.069659 

106.203 

— 106,442 

16 

.892 

,074972 

105,393 

- 106 , 2 ' 5 e 

17 

.956 

.080388 

105,577 

- 105, 960 

18 

1.024 

.086117 

105.133 

- 105,604 

19 

1.096 

.092163 

104,769 

- 105 . 2 o 7 

20 

l.j -72 

.098534 

104,467 

- 105,016 

21 

1.251 

.105232 

104,135 

- 104.352 

22 

1.335 

.112258 

103,213 

- 103,886 

23 

1.422 

i 119614 

102.961 

- 103.607 

24 

l*al 4 

.127295 

102.800 

- 102,913 

25 

1.609 

* 135297 

101,783 

- 102.5 73 

26 

1.708 

.143613 

101.800 

- 102,518 

27 

l.blO 

.152232 

93,633 

- 101,637 

26 

1.916 

,161141 

93,732 

- 100, 771 

29 

2.026 

.170323 

93,330 

- 100.927 

30 

2.138 

.179761 

93.659 

“ 100.719 

31 

2 . 2 s 3 

.189433 

93.?37 

- 100,259 

32 

2.370 

,199314 

96.132 

- 97.857 

33 

2,490 

.209378 

96,576 

- 93,057 

34 

2.611 

.219596 

97,027 

- 98 .A 62 


OacUMEHT SAMPLE CASE 

ThTAjK Tf'TAiRB Ci>PK 

0*000 0*000 0.000000 

O.CJvjO 0*300 0*000000 

j.OoO 0*000 O.OJOOOO 

3.000 0.000 0.000000 

O.OOO 0.000 *BZ‘)95Q 

3.000 0,000 .oZA^tO 

3*030 0.030 . 0 Z 493 Q 

0.000 U.OOO . 82^^350 

3.000 0*000 .BZ 4930 

J.OOO C.OOO 1.234126 


162,119 

- 162.119 

1.234126 

162.204 

- 162.204 

1.234126 

162.267 

- 162.254 

1.234126 

104,235 

- 164.178 

1.234126 

169.304 

- 164,233 

1.526010 

164.460 

- 164,292 

1,526010 

15 7,25 0 

- 157.195 

1.526010 

169.621 

- 164 . 5 C 3 

1.526010 

159,276 

- 159,125 

1*526010 

162,954 

- 162.752 

1.712400 

103.073 

- 162,605 

1.712400 

163.492 

- 155.925 

1.712400 

159,8 15 

- 159.038 

1.712400 

159.839 

- 159.858 

1,712400 

156,369 

- 15^.947 

1.759928 

152,556 

- 156,201 

1 . 75992 b 

149.042 

- 152 , 6 J 9 

1 . 75 S '.28 

l 5 7 . 096 

- 160,480 

1.759928 

160.982 

- 160.380 

1,759928 

ioJ ,9 73 

- 160.400 

1.482451 

157.039 

- 160,504 

1.482451 

161,772 

- 157,207 

1*482451 

161,547 

- 162.337 

1.462451 

157.237 

- 158.161 

1.462451 


CLPK 

AK 

AKOOT 

0.000000 

.080452 

.660795 

O.OuOOOO 

.113164 

.647707 

0 « 0 JOOiiO 

.145226 

.634715 

O.OOOOOO 

■176636 

.621617 

-.OOUOOu 

.207408 

.609011 

-.000000 

.237540 

.596292 

“.OUOOOO 

.26 7039 

.583659 

-.000000 

.295907 

.571110 

— .OOOuOO 

.324151 

.556641 

-aOOOUOO 

.351773 

.546250 

000000 

.378777 

.533935 

-•oooooo 

,405167 

.521693 

-•OOOUOO 

,431277 

■509366 

— (OOOuuO 

.458384 

.496322 

.012951 

•456428 

.482544 

.012931 

.515326 

*463015 

.012951 

■544960 

.452722 

.012951 

.576273 

.436656 

.012931 

.□O 0 O 66 

. 4198 X 2 

.051381 

.637202 

.402167 

,051381 

.068502 

•333703 

. 05 I 3 di 

,699770 

.36 4007 


.031381 .730790 .344671 

•Oslisi . 7 ol 330 .323990 

•081742 .791143 . 3 Q 2 S 87 

.081742 .819973 .260467 

.081742 .647654 .257722 

.081742 .673616 .234327 

.061742 ,697896 .210344 

.082228 .920127 *189616 

.082228 .940059 ,160799 

.082228 .937456 .135340 

*082228 .972106 .109499 

.082223 *983824 .063337 


AKOOT 

-.242741 

-•260766 

-.258662 

-.257035 

-.255243 

-.253504 

-.251616 

-.250160 

-.246592 

-.247052 

-.245559 

-.244112 

-.242691 

-.241221 

-.239725 

-.236191 

-.236 6 30 

-.235049 

-.233457 

-.231661 

-.230272 

-.226699 

-.227152 

-.225641 

-.224161 

-.222763 

-.221454 

-• 2202 M 

-.219056 

-.2iaO0C 

-.217074 

-.216263 

-.215364 

-.215042 



dtfDRe 



N 

XtS) 

Y(N) 

OAitK) 


1 

-.843510 

.508944 

— .J01116 


Z 

-.843511 

.508942 

-.001798 


3 

-.843511 

.508942 

-.002423 



-1,194089 

.600780 

-.j0't710 


5 

-1.Z25271 

,62«499 

-.007114 


6 

-.898120 

,653533 

-.002931 


7 

-.912873 

.499431 

-.0j1593 


a 

-.951854 

,660847 

-.001671 


9 

-1,072383 

,621705 

U041.U8 


10 

-1.193338 

,607523 

-.013316 


11 

-1.208890 

,612132 

-.014769 


iZ 

-1,062167 

.617087 

-.005352 


13 

-.871289 

.518034 

, -.0x0307 


14 

-1.203198 

.624112 

-.019593 


13 

-1.14B552 

.63P894 

-.016383 


16 

-.955113 

.659219 

-.011343 


17 

-.908095 

.675676 

-, 316241 


Id 

-.871560 

.595986 

-.323126 


19 

-.878499 

.563013 

-.021884 


20 

-.956572 

.630661 

-.343648 


21 

-.87«923 

.556931 

-»038d0C 


ZZ 

-.904983 

.515029 

-.344100 


Z3 

-.906310 

.482912 

-.036344 


Z4 

-.900928 

,468802 

-.029424 


25 

-.917299 

.413596 

-.030493 

11 

« 4 

12 - 10 

LI5T1 - 3 

LIST2 -3 

11 

■ 4 

12 > 10 

LIST! • 3 

LIST2 <• 3 

11 

* 4 

12 « 12 

LI5T1 -3 

LIST2 - 3 

11 

■ 5 

12 - 4 

Lion « 3 

LIST2 ■ 3 

11 

- 7 

• 12 ■ 12 

LrSTl « 3 

LI3T2 ■ 3 

11 

■ 8 

j2 ■ 9 

LI9T1 ■ 3 

L1ST2 « 3 


XBIN] 

t8(H) 

SAMS INI 

-.843510 

*508944 

-.001116 

-.843511 

.508942 

-,001798 

-.643511 

.508942 

-.002423 

-1.195089 

,608780 

-.004710 

-1.225271 

•628499 

-.0u7il4 

-.898120 

•653533 

-.002931 

-.912673 

.699431 

-.001593 

-.951854 

.660647 

0016.71 

-I.O 723 B 3 

,621705 

-.004168 

-1,193338 

.607523 

-•0138;16 

-1.208890 

.612132 

-.014769 

-1.082157 

,617067 

-,00«55 2 

-.871269 

.618034 

-.010307 

“1.203198 

,624112 

-,019593 

-1.168552 

.638894 

-*016883 

-.955118 

.659219 

-.011643 

-.908C95 

,6/5676 

-•(#1I241 

-.8 71560 

•595986 

-. 02 D 146 

-.873499 

.563013 

-.021084 

“•95o572 

• 63 O 66 I 

-.0436 >48 

-.676923 

,556931 

-.036800 

-.90<t963 

.515029 

-.044100 

-.9 JoilC 

.4829X2 

-.036344 

-.900923 

•4&6602 

-.029424 

-.917299 

.413596 

-.030493 



AFTtR CliALeSCENCE 


ti 

X{N} 

Y(H) 

GAi1(M 

1 

-.643510 

.508944 

-.J01116 

2 

-.843511 

.508Q4? 

-.001798 

3 

-.843511 

,508942 

-.002423 

4 

-1.204309 

.616660 

-.060002 

5 

-.098120 

.653533 

-*u02931 

6 

-.912873 

.690431 

— .001593 

7 

-.954689 

.659433 

-.012714 

8 

-1.065518 

*618605 

-. 012740 

9 

-.871289 

.610034 

-.010307 

10 

-1,148552 

,638894 

-•016883 

11 

-.900095 

,675676 

-.018241 

12 

-.871550 

,595986 

-.020126 

13 

-.878499 

.563013 

-.021884 

14 

-.956572 

.630661 

-. j43t>48 


-.676923 

.656931 

-.038800 

16 

-.904993 

.515029 

-.044100 

17 

-.906310 

.482912 

036344 

18 

-.900920 

.468802 

-.029424 

19 

-.917299 

• 4135<^6 

-.030493 


XBlHj 

-.843511 
-,d^3511 
-1.2u4309 
-.a^aigo 
-.^I2d73 
-.954689 
-1.06551B 
-.871ZB9 
-1.148552 
-.908095 
—.0 71560 
-.878499 
-.936572 
-.678923 
-.904983 
-.906310 
-.900926 
-.917299 


YBtNJ 

.308944 
.SOo 942 
.508942 
. 6 16 660 
.653533 
.699431 
.699433 
.018605 
.618034 
.638894 
.075676 
.393986 
.363013 
« o 3066 1 
.536931 
.515029 
.<102912 
•468802 
.913596 


GAKBtH} 

—.60x116 

-.001798 

-.002423 

-.060002 

-.002931 

-.001593 

-.012714 

-.012740 

-.010307 

-.010883 

-.U18241 

-.020126 

-.021884 

-.043648 

-.038000 

-.044100 

-.036344 

-.029424 

-.030493 


AFTcR COALcSCtHC£ 


rt 

X(N1 

Y{NJ 

GA.KM) 

XB(K) 

r»(K) 


i 

-.847351 

-.503274 

•OJlilb 

-.643510 

.500944 

-.001116 

2 

-.840460 

-.633618 

.001798 

-.843511 

,508942 

-.0011793 

i 

-.860595 

-.673760 

•002419 

-.84351-1 

.508942 

-•002i423 

4 

-1.306839 

-.460092 

.004729 

-1.204309 

.616660 

— .O 6 C: 0 u 2 

5 

-1.250606 

-.465323 

,007223 

-.898x20 

.653533 

-.002:931 

6 

—s 906181 

-.647009 

,003009 

-.912873 

.699431. 

-.001593 

7 

-.917097 

-.674110 

•OOI 06 O 

-.954689 

.659433 

-.012714 

d 

-.925509 

-.663450 

.001844 

-1.065510 

.618605 

-.012:740 

9 - 

-1.377399 . 

-.459310 

,003909 

-.671239 

.618034 

-.010307 

L(J 

-1,348216 

-,330009 

,012594 

-1. 148562 

.638894 

-.U16I8S3 

11 

-1,464450 

-.132322 

.014066 

-.908095 

,675676 

-.01*241 

l2 

*1.119047 

-.660867 

*006227 

-.871560 

.595986 

-.0201126 

13 

-1,004290 

-.669930 

.010340 

-.878499 

.563013 

-.021^84 

14 

-1,224962 

-.517343 

.017804 

-.956572 

.630661 

-.043648 

1 & 

-.973744 

-.650274 

,0148x3 

-.878923 

.556931 

-.03*1800 

Id 

-.950104 

-.670034 

. 0 IO 556 

-.904933 

.515029 

-.044100 

17 

-,920729 

-.659751 

.017940 

-.906310 

.482912 

-.0361344 

Id 

-.080988 

-.604830 

.030330 

-.906928 

,468802 

-.029424 

19 

-.936216 

-.615499 

,014722 

-.917299 

.413596 

-.036493 

' 

20 

-.892*73 

-.534767 

.020796 




21 

-.940011 

-.553977 

.037453 




22 

-.903225 

-.937937 

.034150 




23 

-.899633 

-.466237 

.025524 




2 't 

-.902786 

-.454924 

.019388 




25 

-.933883 

-.374890 

.0308?! 






z 

CO 

CL 

.026418 

.206230 

0.000000 

•047840 

.300531 

o.ooouoo 

.069859 ^ 

.301502 

.003236 

,098534 

.4?blOO 

.01204? 

,135297 

.439962 

.020436 

.179 761 

.370613 

.020557 

.229937 

. '’1 o24o 

,0445o5 

.282390 

. 137005 

-.002477 

CH 

CY 

CSM 

1.065266 

.058002 

-.157920 


i 


THcTAS 

112.2t.071'5 
I0u.!)&6462 
l(lQtZ07'i<i3 
lO'i. 467^0? 
10i«7()3 39U 
98,658703 
97«437d31 
92,298054 

‘ CSH 
•*«UlCl956 


rHETASI 

-il2,26J7l9 
“100*586482 
“l06*44176U 
“iO5*UX>fe90 
-102,573300 
-iOO, 718952 
-98,464543 
-93.993794 



(namelist input) 


SVTXIN 

KFINAL-40 

ZPERT«.05 

ZPEND-.25 

AATACK-30 

ZLEHGTH-78.272 

ZLND-2.59& 

DEE-7.6 

ZINIT-.Ol 

KVEC-30 

KPRINT-10 

KPRINTP-ZG 

SEND 


CASE 3 


47 



Kours 

KW'''’ 

KCQAL 

4i> 

5 

31 

0 

RLAP1^ 

ftITFOOr 

oou 

PBL 

0. 

0. 

3. 

• 0. 

TPINAL 

OflLT 

D-LT’<T'i 

O^ITHAX 

ijO.OCU 

0,010 

• O'^l 

,1^5 

AATAtK 

srfi^'A 


7LAMDA 

3&.0LK; 

.6 30 

.100 

n.ooo 


PbCtJfi-NT 

Sahple 

CASE 

• 



KRCOAL 

KLOUT 

KPtUNT 

KPPINTP 

iSYH 

LAH 

C 

C 

Id 

20 

U 

0 

PAP 

Pf IX 

PCUAL 

L'GIVC 



0. 

c. 

0. 

X. 



iFlNAL 

71NI7 

ZPER7 

2P6ND 

RC 

XTR 

l.OOC 

.010 

.050 

. 2bO 

.050 

0.000 


cPS 

.100 






it NOTH 


OTAMcrtk 


78.2720 


7.60U0 

ZSTAK 

fiSTAS 

7 

T 

J.UUOO 

.0000 

n.oo'io 

0. 0000 

t.tiOSA 

,5973 

.0200 

,237<J 

d.iSOV 

1,1391 

.0400 

,4757 

^•b'^63 

1.6272 

,0600 

,7135 

b (Zblb 

2.0630 

.0100 

,9514 

7,(>27Z 

2.4479 

,1000 

l.ib92 

9t3V2b 

2.7B30 

,1200 

1.4271 

lO.'JSBl 

3.0693 

.1400 

1.664S 

i2.b23b 

3,3076 

.1600 

1.9026 

1A.U091* 

3.4905 

,1800 

2.140D 

Is.bbAA 

3.6425 

,2000 

2.3764 

1/. 2X9(1 

3,7400 

,2200 

2,5 16 3 

Id. /b!>3 

3,7913 

.2400 

2.6541 

20,3507 

3,(1000 

,2600 

3,0920 

21.9162 

3,8000 

,?800 

3.3296 

23 . lb 

B.tOOO 

.3000 

3.5677 

25.0970 

3.BC00 

.3200 

3,6055 

2o.bl25 

3.POOO 

,3400 

4.0433 

2d,X779 

3 ,11000 

,3600 

4,2312 

29.7934 

a.eooo 

.3900 

4.5190 

3X.303O 

B.flOOO 

.4000 

4,7569 

32.b742 

3.P0Q0 

.4200 

4.9947 

34.4397 

3,8000 

.4400 

5. 2326 

36.0051 

3.1 000 

.4600 

5,<r?04 

37.5706 

3,8000 

.4900- 

5.7063 

39.1360 

3,8000 

.OOOO 

■5.9461 

40.7014 

3.8000 

.szon 

6. lb 39 

42.2b69 

3.8000 

.5400 

6,421b 

43,6323 

3,8000 

,5600 

6.o59o 

45,397b 

3.8000 

.5 800 

6.6975 


DUCUMfNF CASc 

dlDY GiO'ltTkY 



ARt A 


8 


AVi 


45*3646 


10.2/09 


3*8000 

RZkQ 

OrcZRti 

DClRZKO 

4 

adot 

AODOT 

• OJOU 

6.23J3 

-37.3531 

.OOuO 

• 6923 

-.2677 

.1572 

7.4900 

-36.5125 

.1572 

,6298 

-.2562 

,2998 

6. 7 716 

-35,3507 

• 2vv3 

,5 694 

-.2500 

.4232 

6.0749 

-34.3460 

,42t2 

.5 108 

-.2429 

.5 429 

5,3969 

-33,4607 

»5429 

■ 453b 

—•2367 

.u4(t2 

4.734b 

-32. 7406 

.6442 

• 3961 

-.2315 

.7324 

4.0365 

-32,1136 

.7329 

. 3436 

-.2271 

,o077 

3.4496 

-31.5909 

.8077 

.2901 

-,2234 

,67o4 

2.U22Z 

-31.1644 

,8 704 

.2373 

-.2204 

.9207 

2,2024 

-30.6262 

,92o7 

.185 2 

-,2iev 

.9536 

1,5665 

-30,5777 

«95oo 

.1336 

-.2162 

.9642 

.973 8 

-30,4097 

,9842 

.0823 

-,216w 

.9977 

.3716 

-30,3217 

;9977 

.(/312 

-,2l44 

1.0000 

C.uOOC 

0.0000 

l.OOJd 

O.JOOO 

0,0000 

1,0000 

0, OUOLi 

O.OOUC 

l.OOOO 

0*000(J 

0,0000 

1.0C30 

0,0000 

0.3000 

l.OwoO 

0.0000 

0.0000 

1.0000 

0,0000 

0.0000 

i^OOCO 

o.oouo 

u.oooo 

1.0000 

0, UOOO 

O.OOuO 

l.Oooo 

0,0000 

o,dOO(< 

l.OOOO 

0,0000 

o.OOuO 

1.0000 

0,0000 

OaOUbO 

i.uijOO 

0 , oooc 

0.0000 

1 tOiioo 

u.oooo. 

0,0000 

1.0000 

0,0000 

0.0000 

1.0000 

o.oouo 

0,0000 

1.0000 

0.0000 

O.OOoO 

1.0000 

O.DOCO 

O.OuCO 

l.UoOO 

0. OOOO 

0.0000 

1,0000 

o.oouo 

o.ooou 

1.0000 

o.oooc 

o.oooc 

1,0000 

o.oouo 

o.ooou 

1.0000 

0,0000 

0. JOOC 

1,0000 

O.uOOU 

' o.ooco 

1 .uOJO 

o.oooc 

0,0000 

1,0000 

0. 0000 

J.UuOO 

l.OOuO 

(l.OOOO 

• 0.0 OuC 

1,0000 

J.UuUO 

0, OObu 

l.OOOO 

o.ooou 

o.oooc 

1, OOuO 

o.uooo 

o.ooou 

1.0000 

0.0000 

O.OOUO 

ifCOuO 

U.OOUu 

O.UOOO 

1 .0.1 JU 

O.oooc 

O.OOOO 

l.OOuO 

0,0 OQO 

O.UOOu 


oi 

o 


■^6.9632 

3.8000 

,6000 


3.ROOO 

.6200 

&0.U9^i, 

3.8000 

,6400 

51.6596 

3.8000 

.6600 

53.2250 

3.8000 

-.60 00 

5<r,79u4 

3.8000 

.7001 

56.3556 

3,8000 

.7200 

57,9213 

3,8000 

.7400 

5v.4o67 

3.8000 

,7600 

61,0522 

3.8000 

.7800 

62.0176 

3.0000 

,8000 

64.1830 

3.8000 

.8200 

66,7485 

3.8000 

.9400 

6 7,3139 

3.»-;ouo 

.8600 

60, 0794 

3.POOO 

. .8800 

70. <,440 

3.8000 

.9000 

7<i,Ul02 

3.8003 

.9200 

73,6757 

3**1600 

.9400 

75.1411 

3.8000 

.9600 

76.7066 

3.-J000 

.9000 

78.^720 

3. 8 000 

l.OOOo 


7,1353 

I.OOCO 

0,0060 

7.3732 

1.0000 

0.0000 

7,bllO 

1.0000 

0. 0000 

7.8469 

l.OwUC 

O.OOOt 

6 * 656 7 

l.OOUO 

0.0000 

H.3245 

l.oOOO 

0,0000 

8.5624 

i.ouoc 

0.0000 

B.O0U2 

l.UCrOO 

O.OOob 

9.6381 

I.OUOO 

O.oOOO 

9.2759 

l.UUOO 

O.OOOO 

9.513b 

l.UOUO 

0.0000 

9, 7516 

l.OO'JO 

O.OOOO 

9.9894 

1,0000 

0.0000 

10.2273 

1.0000 

o.ooao 

10,4651 

1.U3O0 

0. JOUC 

10* 7i]30 

l.UODO 

O.OuOO 

10,9468 

1. JooO 

0. OuOO 

11 . 17t7 

l.uk/’Jb 

o.uOOO 

11.4166 

1.0000 

O.OCOu 

11.6544 

l.uOOO 

O.OOUU 

11.8922 

l.lMiOO 

O.oOOO 


O.OOOO 

l.uuu 0 

b.uQOO 

0.000b 

1.0003 

O.OOOO 

0. oOOO 

,1.0000 

o.cooo 

0.0000 

T.ouoo 

O.oOOO 

U.-JOOO 

1 . 0 JuO 

O.uOOU 

0.0000 

I.OUOO 

0.0000 

0,0000 

l.OOoO 

O.OOOO 

O.OOOO 

i.Oooo 

0.0000 

c.oooc 

l^OOvfO 

O.OCOU 

0,0000 

l.UOUO 

u.uooo 

O.OoOO 

1 .OOo J 

O.OOOO 

0.0000 

1 .0 Jo 3 

O.OOOO 

0.0000 

IfOOJO 

O.UOOu 

0. 3000 

l.OOOO 

0. OuUU 

0.0000 

1.0 buO 

O.OOOO 

0.0000 

l.uooo 

0 . w OuO 

o.oujo 

1*0000 

U.OOOu 

0. OObC 

l.COO J 

0.0000 

O.OCOU 

l.OuOl) 

u* OuOb 

0. OubO 

IfOooo 

0 . oOOu 

O.OOOO 

1*0000 

w* OObo 


tiadUbU 

OfOCfOu 

0 «UuUU 

0.0000 

O.OQOO 

0.0000 

C.oouo 

o.ooou 

O.QOUU 

0.0000 

O.OObU 

j.aoou 

O.OUUO 

0.0000 

O.OOCiu 

O.oooo 

O.UllOu 

0 * 0 o 0 u 

0.0000 

O.OuOO 


CONTINUED FROM PREVIOUS PAGE 



DOCUK£NT SAfiPlt CkSi 


11 Ht 

Z 

TH=TA9 

THPTASB 

THTASR 

THTASRb 

COPK 

CLPK 

AK 

AKDCiT 

AKOOT 

■ 0.19 

■OlOCOO 

109.000 

-130,000 

O.OuO 

0 . QOv 

0.300000 

O.OQUOOO 

•U60952 

,660795 

-.262798 

• 169 

.019209 

199.399 

-199.399 

0.000 

u. 000 

0.000000 

O.OUOwOO 

.U3169 

.697707 

-.260786 

■ 219 

,016909 

117.917 

-117.917 

0 • ooo 

0. DuO 

0.0^2^000 

0.000000 

.195225 

•039715 

-.258682 

• 269 

.022613 

119.720 

-119. Uh 

O.OuO 

C.OOO 

O.OOoOCC 

0.000000 

• 176630 

.621817 

-.257035 

.319 

,0269ia 

112.261 

-112.261 

0.000 

0. 000 

.P29950 

uouooo 

.20Z906 

.609011 

-.255293 

■ 369 

.031022 

111.530 

-Ut.630 

o.uuo 

0.000 

.629950 

-•OJOOwQ 

,237590 

.596292 

-.253509 

.919 

.035227 

110.736 

-110. 7b6 

0 . OOO 

0.000 

.629950 

-.000000 

•267039 

.583659 

-.251616 

*969 

.039931 

110.197 

-no. 197 

O.uoO 

0.000 

.629950 

-•(jQOOOu 

.295907 

.571110 

-.250160 

.619 

.093635 

109,259 

-109.259 

0.000 

0.000 

.629950 

-.000000 

.329151 

.558691 

-.291592 
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OUCUMtNT iArtf'Lt GASU 


K 

Tlrtu 

1 

TH=TA9 

THrTASb 

UfAiR 

ThTASRB 

CORK 

tLRK 

AK 

AKD01 

AKOUT 

iCi 

• 669 

.O97«A0 

10=^,5*»o 

“193 , atj6 

o.ooj 

0.000 

1,239126 

— .00 OOOli 

,351773 

.596250 

-.297052 

11 

• Qi9 

.052099 

lOR . l?3 

-133,123 

lo2.xl9 

-162, 119 

1.239120 

-.OOOOOw 

.376777 

.533935 

-.295559 

12 

• 66v 

.U5f299 

107.965 

“137.9tob 

,62 • 209 

-162.209 

1.239126 

-.uOOOuO 

,905167 

.321693 

-.299112 

13 

.720 

.0605 OH 

106.967 

-107,033 

162.267 

-162.239 

1.239x26 

-.UOtOoO 

.93X277 

.509360 

-.292691 

1^ 

• 77i» 

.066C91 

106.959 

-136.620 

109.235 

-169,178 

L. 239x26 

— • UOOOUO 

. ‘*56389 

.996322 

-.291227 

lb 

• o3i 

.069^69 

106.203 

-106,692 

169.309 

-169. 233 

1.526J10 

,012951 

. ■f6692U 

.962599 

-.239725 

16 

• U92 

.079972 

105,993 

-106. 2bo 

169.900 

-169.292 

1.5260X0 

.012951 

• 5x5326 

•9630x5 

-.238191 

17 

.956 

.080329 

105.577 

-105,960 

167.250 

-157.195 

1.526010 

.012951 

.599980 

.952722 

-.236630 

id 

i.u29 

.086117 

195.183 

-105,609 

16^.621 

-169.303 

1.5260X0 

,012951 

.575273 

.936656 

-,235099 

19 

1.096 

.092163 

109.769 

-105.267 

15 ) . 2 76 

-169,125 

1.526010 

,012951 

*d0o066 

.9198X2 

-.233957 


H 

K 

X 

y 

xonx 
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5.0 OPERATION 


The following example shows a deck setup for the execution of 
DIV0RCE. The cards designated by a plus sign (+) are those which are 
optional depending upon the existence of special inputs. The job 
and any accounting cards are not shown. The setup shown assumes 
the total use of disk files. If tapes are to be used the proper 
control cards must be substituted. Also the input data is assumed to 
be on cards rather than a disk file. 


A 


ATTACH, DIV0RCE. 
+REQUEST, TAPES ,*PF. 
+REQUEST, TAPES, *PF. 
+REQUEST,TAPE7,*PF. 
+ATTACH ,TAPE4 , RESTART . 
ATTACH, RDALIB. 
LIBRARY,RDALIB. 
LDSET,PRESET=7.ER0,MAP= 
DIV0RCE. 

+CATA10G, TAPES, NAME. 
+CATAL0G, TAPES, NAME. 
+CATAL0G,TAPE7,NAME, 
E0S. 

TITLE CARD. 

$ ENGAGE 


Binary version of the program 
If RITE0UT=T 
If KL0UT?«O 
If KPRESjFSO 
If READIN=1 

Vector routine library 

Set core to zero and map option 
Execute. 


End of section 
Case title card 
NAMELIST title card 
Data 


$END End of case card 

E0I End of infomiation 

If KVEC is large enough that the vector routines are not being used, 
the RDALIB cards may be left out, although load errors (unsatisifed 
externals) will be generated. These are generally non-fatal. Core should 
be set ot zero for safety since the program was written on a system that 
did this automatically. 
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